function y = func_alfa_2(p, CONSTS)
    eps = CONSTS.eps;
    g = CONSTS.g;
    q2 = func_q2_from_p(p, CONSTS);
    
    y = -1 + (p.^2 + q2.^2 - eps)/g;
end
